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Abstract. I present results from a set of 3D spherical-shell MHD simulations 
of the buoyant rise of active region flux tubes in the solar interior which put 
new constraints on the initial twist of the subsurface tubes in order for them to 
emerge with tilt angles consistent with the observed Joy's law for the mean tilt 
of solar active regions. Due to the asymmetric stretching of the fi-shaped tube 
by the Coriolis force, a field strength asymmetry develops with the leading side 
having a greater field strength and thus being more cohesive compared to the 
following side. Furthermore, the magnetic flux in the leading leg shows more 
coherent values of local twist a = J • B/£? 2 , whereas the values in the following 
leg show large fluctuations and are of mixed signs. 



1. Introduction 

If we believe that active regions on the solar surface originate from a strong, 
predominantly toroidal magnetic magnetic field generated at the base of the 
solar convection zone, then we need to explain how the magnetic flux tubes rise 
through the solar convection zone to the surface. Despite the turbulent nature 
of solar convection, bipolar active regions on the solar surface show remarkable 
order and organization as described by the Hale polarity law and Joy's law 
of active region tilts. This suggests that the field strength of the rising flux 
tubes are significantly super-equipartition compared to the kinetic energy of 
solar convection in the bulk of the solar zone such that emerging flux tubes are 
not severely distorted by the convective motion and roughly preserve the original 
orientation of the toroidal magnetic fields at the base of the solar convection 
zone. 

Previously a large body of calculations based on a highly simplified thin 



flux tube model (|Spruitlll981l ) has provided important insights into the dynamic 



evolution of a Q-shaped rising flux tube in the solar convective envelope (e.g 

D'Silya and Choudhuri Hl993l;lFan. Fisher, and DeLuca II 19931: Caligari. Moreno-Insertis. and Schiissler 



19951 ; ICaligari. Schussler. and Moreno- Insertis 1998l: lFan and Fisher 1119961 ). These 



calculations suggest that the field strength of the toroidal magnetic field at the 
base of the solar convection zone is in the range of about 3 x 10 4 G to 10 5 G 
in order for the emerging tubes to be consistent with the observed properties 
of solar active regions. It is found that the Coriolis force acting on the ris- 
ing ^-shaped tube can produce asymmetries between the leading and following 
legs of the fi-tube which provide explanations for several observed asymmetric 
properties between the leading and following polarities of bipolar active regions. 
These include (1) a slight tilt of the emerging loop consistent with the mean 
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Figure 1. The evolution at the tube apex of (top-left panel) the Alfven 
velocity va, ris e velocity iy, convective vel ocity v conv based on JCD's solar 
interior model (jChristensen-Dalsg aard 1982), the azimuthal velocity v$, (top- 
right panel) the magnetic field strength B, (bottom-left panel) time elapsed 
since the onset of the Parker instability, and (bottom- right panel) the tube ra- 
dius a and the local pressure scale height H p , as a function of depth, resulting 
from a thin flux tube simulation of an emerging ^-shaped tube described in 
iFan and Gong I (|1999t corresponding to the case shown in the top panel of 
Figure 1 in that paper) 

tilt of bipolar active regions; (2) the asymmetric geometric shape of the emerg- 
ing loop which may give rise to the observed asymmetric horizontal motions of 
the two polarities of an emerging active region, and (3) the asymmetry in the 
field strength between the two sides of the loop with a stronger field along the 
leading side which may be the cause of the observed asymmetric morphology for 
the leading and following polarities. Since the thin flux tube model satisfies the 
frozen-in condition perfectly, it faithfully describes the evolution of the magnetic 
field strength, the magnetic buoyancy, and hence the velocity of the rising flux 
tube without being subject to numerical diffusion. Figure 1 shows the evolution 
of various quantities at the apex of the rising tube as it traverse the convection 
zone, resulting from a thin flux tube simulation of an emerging O-shaped tube 
developed due to the Parker instability of a 10 5 G toroidal flux tube initially 
in mechanical equi librium at the base o f the solar convection zone at an initial 
latitude of 15° (see lFan and Gong~lll999l . for details of the simulation). It can be 
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seen that the rise velocity v r remains quite small (< 200ms" 1 ) and the Alfven 
speed v a remains nearly constant (being much greater than both the rise and the 
convective flow speed) in the bulk of the convection zone, until the top few tens of 
Mm of the convection zone, where v r accelerates steeply and v a decreases rapidly 
due to the steep super-adiabaticity in this top layer. At a depth of roughly 20 
Mm, the radius of the tube a exceeds the local pressure scale height H p and v r 
also exceeds v a . At this point, the thin flux tube approximation breaks down 
and the tube is likely to be severely distorted and fragmented. Nevertheless, if 
one continues to use the v r from this point on as an estimate, one finds that the 
tube will rise through the last 20 Mm depth of the convection zone in a time 
of about only 7 hours. These results provide some basic information for local 
helioseismology to estimate possible helioseismic signatures (e.g. acoustic travel 

time perturbations) for detectin g subsurface emerging active region flux tubes. 

Later MHD simulations (e.g Emonet and Moreno-Insertis Ill998l ; Abbett. Fisher, and Fan 

20001 1 of buoyantly rising flux tubes in local Cartesian geometries suggested that 



in order for the flux tube to rise cohesively, a minimum twist of the flux tub e 
field lines given by q = BQ/(rB z ) = (H p a)~ 1 / 2 is needed ( Longcope et al.|[l999l ). 
where q corresponds to the angular rate of field line rotation about the axis 
over a unit length along the tube, B z and Bq denote respectively the axial and 
azimuthal magnetic field in the cylindrical flux tube, r is the distance to the 
tube axis, H p is the local pressure scale height, and a is the tube radius. This 
minimum twist is found to be about an order of magnitude too big compared 
to the twist deduced from yector magnetic field observations of sol ar ac t ive re - 
gions on the photosphere ( Longcope et al.lll999l ). More recently, iFan I (2008) 
has carried out a set of 3D anelastic MHD simulations of the buoyant rise of 
active region scale flux tubes in the solar interior in a rotating spherical shell 
geometry. These simulations put further constraints on the twist of the interior 
flux tubes in order for them to emerge with tilts that are consistent with the 
observed mean tilt ang les of solar active regions. Cons istent with previous thin 
flux tube calculations ( Fan. Fisher, and DeLuca 19931 ). it is found that due to 



the asymmetric stretching of the flux tube by the Coriolis force, the O-shaped 
emerging tube develops a field strength asymmetry with the leading side having 
a stronger field, and thus being less fragmented compared to the following side. 
We will review these results in the following and further examine the asymme- 
try in the field line twist and magnetic helicity transport along the leading and 
following side of the emerging tube. 



2. Results from 3D anelastic spherical-shell simulations: 



Fan ( 20081 ) performed a set of 3D anelastic MHD simulations of the buoyant rise 



of magnetic flux tubes in a rotating spherical shell domain given by: r = [r c , rt] 
where r c = 0.722i? Q is the base of the convection zone and r t = 0.977-R© is at 
about 16Mm below the photosphere, 9 = [0,7r/2], and <fi = [0,7r/2]. The domain 
is resolved by a grid with 256 grid points in r, 512 grid points in 0, and 512 
grid point in (j). The grid is uniform in r and 0, but nonuniform in 8, with finer 
grid spacing in the region traversed b y the rising flux tube and co arser grid in 
the polar region. JCD's solar model ( Christensen-Dalsgaardl Il982l ) is used for 



the reference stratification in the simulation domain, and the superadiabaticity 
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in the convection zone is set to zero, so that the fluid is marginally stable to 
convection. Thus we are modeling the dynamic rise of buoyant flux tubes in a 
quiescent model convective envelope without the effects of the convective flows. 
Initially a twisted toroidal flux tube with an axial field strength of 10 5 G and 
a radius of a = 5.59 x 10 8 cm ~ 0.1H P is placed just slightly above the base of 
the model convective envelope. For flux tubes with such a strong initial field 
strength, which is 10 times the field that is in equipartition with convection, 
the dynamic effe cts of convective flows on the buoyant flux tube is expected 
to be small (e.g. iFan. Abbett. and Fisher" 20031 ) and neglecting convection in 
the simulations is expected to be a good approximation. An initial sinusoidal 
variation of entropy is imposed along the toroidal tube such that the middle of 
the toroidal tube is approximately in thermal equilibrium with the surrounding 
and thus is buoyant and the two ends (at <f> = and <j) = 7r/2) of the toroidal 
tube are in neutral buoyancy. Thus an $7-shaped rising flux tube develops due 
to the initial buoyancy specification. 

Simulations where the initial twist of the field lines is set to the critical value 
needed for a cohesive rise of the tube, q = ±0.3a- 1 « ±(H p a)- 1 / 2 , show that 
the final tilt of the fi-tube at the apex is largely determined by the sign of the 
twist, with positive (negative) q resulting in a clockwise (counter-clockwise) tilt. 
Here q denotes the angle of field line rotation about the axis over a unit distance 
along the tube. In these cases, it is found that the writhing of the tube due to 
the field line twist dominates the effect of the Coriolis force and produces the 
wrong sign of tilt at the apex of the emerging tube, if the observed hemispheric 
preference of the sign of twist (i.e. left-handed-twist or negative q in the northern 
hemisphere) is assumed. It is found that in order for the emerging tube to show 
the correct tilt direction (consistent with observations), the initial twist of the 
flux tube needs to be less than a half of that needed for a cohesive rise. Under 
such conditions, severe flux loss is found during the rise. Figure 2(a) shows the 
3D volume rendering of the magnetic field strength of an O-shaped rising flux 
tube whose apex is approaching the top boundary, resulting from a simulation 
with an initial twist q = — 0.15a -1 , initial field strength of 10 5 G, starting from 
an initial latitude of 15°. Only about 45% of the initial flux remains in the 
coherent portion of the tube that rises to the upper boundary. In this case the 
final fi-tube attains a clockwise tilt angle of about 6.5° for its apex portion (as 
can be seen in the cross-section shown in Figure 2(b)), consistent with the mean 
tilt of solar active regions. Furthermore, the Coriolis force drives a retrograde 
flow along the apex portion, resulting in a relatively greater stretching of the 
field lines and hence stronger field strength in the leading leg of the tube. With 
a greater field strength, the leading leg is more buoyant with a greater rise 
velocity, and remains more cohesive compared to the following leg. Figure 2(c) 
shows selected field lines threading through the coherent apex cross-section of 
the rising J7-tube. One can clearly see the asymmetry where the field lines in the 
leading side are winding about each other smoothly in a coherent fashion, while 
the field lines in the following side are fraying out. We have also evaluated the 
local twist rate given by a = J • B/£? 2 , where J is the electric current density, 
along each of the selected field lines shown in Figure 2(c) and plotted them (as 
dots) as a function of depth in Figure 3 for the following side (left panel) and 
the leading side (right panel) respectively. The mean value averaged over the 
field lines as a function of depth is shown as the solid curves. We can see that 
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Figure 2. (a) 3D volume rendering of the magnetic field strength of an 
rising f2-tube whose apex is approaching the top boundary, resulting from a 
simulation with an initial twist q = —0.1 5a -1 , initial field strength 10 5 G, and 
an initial latitude of 15° (see the LNT run in iFan N2008I ): (b) A cross section 
of B near the top boundary at r — O.937i?0; (c) selected field lines threading 
through the coherent apex cross-section of the fi-tube. 
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Figure 3. Dots show values of a = J • B/£? 2 computed along each of the 
selected field lines of the final Omega-tube shown in Figure 2(c) as a function 
of depth for the following side (left panel) and the leading side (right panel). 
The field line averaged mean a is shown as the solid curve. 
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the a values along the field lines on the following side show large fluctuations 
and have mixed signs, whereas the values along the field lines on the leading side 
are more coherent and are consistently negative. However the mean a averaged 
over the field lines does not show a significant systematic difference between the 
leading and following sides. 



3. Conclusions 

Our 3D simulations of buoyantly rising twisted flux tubes in a rotating spherical 
shell geometry suggest a new upper limit for the twist of the subsurface flux 
tubes in order for them to emerge with tilt angles consistent with the observed 
tilts of the majority of solar active regions. This upper limit is about a half of 
the twist needed for a flux tube with the magnetic buoyancy to rise cohesively. 
The simulation results also show that due to the Coriolis force, the leading 
leg of the fi-shaped tube tends to have a stronger field strength, being more 
buoyant and more cohesive compared to the following leg. Furthermore the field 
lines in the leading leg show more coherent values of local twist as measured by 
a = 3 ■ B/5 2 , whereas the field lines in the following leg are more frayed and 
sh ow larger fluctuations a nd m ixed signs of twist. Recent observational studies 
by iTian and Alexander ( 20081 ) have found an systematically greater helicity 



injection flux in the leading polarity than the following in several emerging active 
regions. Although in our model, the field lines of the leading leg of the O-tube 
do not show a systematically greater twist on average, the greater buoyancy 
and hence higher rise velocity of the leading leg in conjunction with its more 
coherent twist may give rise to a greater upward helicity flux in the leading 
polarity comparing to the following as a result of the emergence of the O-tube. 
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